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ABSTRACT 

We investigate the early-time light-curves of a large sample of 223 type II supernovae 
(SNe) from the Sloan Digital Sky Survey and the Supernova Legacy Survey. Having 
a cadence of a few days and sufficient non-detections prior to explosion, we constrain 
rise-times, i.e. the durations from estimated first to maximum light, as a function of ef¬ 
fective wavelength. At restframe g'-band (AeS = 4722A), we find a distribution of fast 
rise-times with median of (7.5 ±0.3) days. Comparing these durat ions with analytical 
shock models ofiRabina k &: W axmanI (l2013l l: iNakar fc Saril (l20I0ll and hydrodynami- 
cal models of iTominaga et al.l (l2009ll ~ which are mostly sensitive to progenitor radius 
at these epochs, we find a median characteristic radius of less than 400 solar radii. The 
inferred radii are on average much smaller than the radii obtained for observed red 
supergiants (RSG). Investigating the post-maximum slopes as a function of effective 
wavelength in the light of theoretical models, we find that massive hydrogen envelopes 
are still needed to explain the plateaus of SNe II. We therefore argue that the SN II 
rise-times we observe are either a) the shock cooling resulting from the core collapse of 
RSG with small and dense envelopes, or b) the delayed and prolonged shock breakout 
of the collapse of a RSG with an extended atmosphere or embedded within pre-SN 
circumstellar material. 

Key words: supernovae - red supergiants. 


1 INTRODUCTION 

In the standard picture of single massive star evolution, red 
supergiant stars (RSG) with zero-age main-sequence mass 
between approximately 8 and 25 Mq end their lives with the 
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collapse of their cores (CC) in luminous explosions known 
as ty pe II supernovae or SNe II (for a review see e.g. ISmartd 
I 2 OO 9 II . The progenitors retain part of their hydrogen layers 
prior to explosion, so that their SNe display hydrogen-rich 
spectra and varying light-curve shape: from linearly declin¬ 
ing after maximum, commonly known as SNe IIL, to plateau 
supernovae, SNe IIP, with characteristic post-maximum 
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plateau dBarbon et alJ 1 19791 : iHeeer et alJ l2003ll . The vari¬ 
ation of plateau lengths is thought to be mainly produced 
by the amount and density of hydrogen which produces the 
reco mbination wave traveling inwards in mass coordinates 
Litvinova fc NadezhinI Il985l : iKasen fc Woosle^ l2009l : 
Bersten et al.ll2011f l. Recently, the study of large samples of 
SNe IIP and SNe IIL have driven discussion as to whether 
these two groups a c tually form a single con t inuous set (e. g. 
Arcavi et alJ [ 2 OI 2 I: lAnderson et alJ l2014l: iGutierrez et alJ 


2008; ISchawinski et af] l2008l : SNLS-06Dljd iGezari et alJ 


2 OO 81 I . The detection of this flash can give tight constraints 


2 OI 4 I : ~anders et alJ 2014 : Faran et al. 1201^ al. In addi¬ 

tion, two other hydrogen-rich SNe are identified: SNe Iln 
and SNe lib. SNe Iln are supernovae whose ejecta inter¬ 
act with dense circumstellar material producing distinctive 
narrow hydrogen emissio n lines in their spectra and pow¬ 
ering the ligh t-curve fe.g. ISchlegellll99nl : iTaddia et al.ll201,ll : 
iMoriva et al.ll201lh . SNe Ilb belong to the class of stripped- 
envelope explosions. They first show hydrogen-rich spectra 
similar to SNe IIP b ut then evolve to hydrogen-deficient 
SNe l b at later times llFilippenko et al.lll993l : iNomoto et al.l 

ll99,ll : [^ufano et al.ll2014h . In this paper we focus on SNe IIP 

and SNe IIL, calling them together “SNe 11”. 

The best evidence for SNe II originating from RSG 
stem s from their direct id e ntification in pre-explosion images 
(e.g. ISmartt et all l2009l : iFraser et al.l I 2 OI 2 I : iMaund et al.l 
l2013h. although this was predicted by early theoretical works 
dGrassberg et al.lll97ll : IChevalierll 19761 : iFalk fc Arnett|[l977l : 
lArnettlll980l L These studies possibly indicate a mass range 
for SN II explosions of 8-I8M0, which may suggest that 
more massive (observed) RSGs do not explode as lumi¬ 
nous events. Thi s has been referred to as the “red super¬ 
giant problem” dSmartt et akllioogl L Low upper limits on 
the initial mass have been inferre d from nucleosynthetic 
yields from nebular spectra analysis djerkstrand et aDl20l4 
iDessart et akllioi^. however hydrodynamical models yield 
large r masses ( Utrobin fc Chugall20ok l2009l : iBersten et al.l 
I 2 OIII L On th e other hand, the fa mily of type II S N 1987A- 
like SNe (e.g. iTaddia et al.ll201^ : IPastorello et al.ir2012h are 
well explained with b l ue supergiants (BSG) progenitors 
dShigevama et al.l 1 19881 : IWooslev et al.l 1 19881 : lArnett et al.l 
1989|). Compared wi th radii of RSG (200 < R/Rq < 1500, 


Levesc[ue_jt_^ 200^ , BSG are much smaller [R/Rq < 200, 
and are in principle not expected to be 


Kudritzki et al. 


200i 


end points of single stellar evolution. However the specifics 
of e.g. convection, mixing, magnet ism and binar y interaction 
can allow for such explosions (see lLangerj[2012l L 

A key ingredient to study the properties of the super¬ 
giants that result in SN II explosions is the early light-curve 
prior to the characteristic plateau where hydrogen recom¬ 
bination dominates the SN display. The first electromag¬ 
netic radiation emitted by a SN occurs when the shock 
wave created from CC emerges from the stellar surface (e. g . 
Colga^ 197^ Falk fc Arnettl[l97^ : iKlein fc Chevalie'rlll97d : 


Ensman fc Burrows! 1992l L This shock breakout emission 


is expected to be of short duration (few hours for typical 
SNe IIP) and has been the subject of recent studies as new 
observing capabilities offer the possibility to detect it. It is 
brightest in high frequencies. X-ray and ultraviolet (UV), 
and a few claims of det ections exist for some hydrogen- 
free CC SNe (SN 200 6aj, ICampana et al.l [20061 : SN 2008D 
ISoderberg et al.l l2008l L and for two SNe IIP with optical 
data from the Supernova Legacy Survey and near-UV from 
the Galaxy Evolution Explorer fSNLS-04D2dc. rGezari et al.l 


on progenitor char acteristics (e. g. iTominaga et al.l l2009l : 
iNakar fc Saril l201Cll L The posterior rise is powered by the 
shock cooling as the heated expanding stellar envelop e dif¬ 
fuses out (e. g. iGrassberg et aLlll97ll : rChevalietlll976ll . The 
study of the early rise of the main light -curve can also giv e 
revealing constraints on the progenitor (iGezari et al.ll2C)l(]l L 
Despite the growing number of nearby SNe II discov¬ 
ered with high quality data, the early parts of the light- 
curves are still quite unexplored, yet fundamental to con¬ 
strain progenitor characteristics. Some constraints on in¬ 
dividual well observed nearby SNe II exist and show that 
they rise faster than 1 0 days in B-band and more slowly 
at longer wavelengths |Leonard et al.|j2002al f9: ISahu et al.l 


l2006l : [Tomasella et al.|[2013l L Recentlv. lGall et a,iT ( 2015lV an- 

alyze the rise-times of 19 low-z SNe II conhrming short rise- 
times and Ending some evidence for a relation between the 
rise-time and the peak brightness. 

In this work we study the early light-curves of a large 
sample of SNe II (SNe IIP and SNe IIL) from two rolling 
searches with good cadence and sufficient data prior to ex¬ 
plosion: the Sloan Digit al Sky Survey-II Supernova Sur¬ 
vey (hereafter SDSS-SN. ISako et al.ll2014ll and the Super¬ 
nova Legacy Su rvey (hereafter SNLS, Astier et al.l I2OO6I : 
IGuv et ^ l20ld j. By comparing the data with analytical 
shock models and hydrodynamical shock models, we are able 
to constrain specihc progenitor properties. A plan of the pa¬ 
per follows: in section [5] we shortly review the theoretical 
models used in the analysis, in sectionOwe present the light- 
curve data, in section|4]the tools we use in the analysis, in[5] 
we show our main results and summarize them in[^ 


2 SHOCK MODELS 

2.1 Analytical shock models 

Analytical models for the shock breakout and immediate 
evolution of an exploding supergiant star offer the advan¬ 
tage of being more practical than more complicated nu¬ 
merical models for htting light-curves with varying pro¬ 
genitor parameters in an unexpensive way. However their 
treatment is simplistic and need to be tested with more so¬ 
phisticated calcu lations (see se c tion JpJ. We use the ana- 
lytical models o fiNakar fc Saril (l201ol . hereafter NSIO) and 
iRabinak fc WaxmanI ( 2011 . hereafter RWll). These models 
provide analytic expressions relating the emission and the 
model parameters such as progenitor mass, radius, and ex¬ 
plosion energy. Both models assume pre-explosion density 
prohles of p{r) oc (R* — r)", where R* is the progenitor ra¬ 
dius and the power-law is n = 3 for radiative and n = 3/2 
for convective envelopes, typical of BSG and RSG stars re¬ 
spectively. The assumed density prohles are only valid for 
the outermost part of the progenitor structure. The models 
start post explosion, with that of NSIO including the actual 
shock breakout, and follow the spherical expansion phase of 
the stellar envelope until the onset of recombination and/or 
radioactive decay. They can thus be used in principle, de¬ 
pending on the model parameters, for the hrst hours up to 
several days for red supergiants (see eq. 16 in RWll). After 
this, when recombination sets in, the assumption of con¬ 
stant opacity starts to break down. RWll include a model 
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Figure 1. Left: Early photospheric temperature and bolometric luminosity evolution for models of NSIO and RWll of a red supergiant 
progenitor with 500-R©, ISM© and 1 foe (10®^erg). The epoch is with respect to explosion. Right: Fluxes in filters NUVgalbx and 
Sloan griz for the same models. Stars show three characteristics times (see text). 


extension for varying opacities but we do not use it in this 
analysis. 

The analytical models of RWll and NSIO can be seen 
in Figure [T] for a standard scenario with R = 500/?©, 
M = ISM© and E = 10®^erg (1 foe). The bolometric lumi¬ 
nosity and photospheric temperature evolution of the first 10 
days are shown, as well as the corresponding prediction for 
NUV and optical bands. The NSIO models show in all bands 
a radiation precursor or shock breakout, i.e. an initial peak 
and decrease at very early phases, followed by the usual rise 
to maximum light. This behaviour originates from the evo¬ 
lution of the spectral energy distribution as the bolometric 
luminosity drops together with the temperature as well, so 
that the tail of the radiation at longer wavelengths is higher 
at later times due to lower temperatures (see Figure [2)). The 
duration and shape of the precursor as well as the rise of 
the light-curve depend on model parameters, particularly 
the radius, as will be shown later. RWll do not model the 
very early phases of shock breakout. 


2.2 Hydrodynamical models 


Numerical calculations to provide models for the early evo¬ 
lution of the emission of a SN have been developed by sev- 


eral authors dFalk & Arnettlll977l:lEnsman & Burrows! 

1992; 

Blinnikov et al.l 2000l: Utrobin 

2 OO 7 I: Gezari et al. 

2008; 

Tominaga et al.ll2009l: iBersten et al.ll201l|). Thev have the 


advantage with respect to analytical models of treating the 
pre-SN stellar structure and the physics in a more realistic 
mann e r. We focus here on the models by iTominaga et al.l 
ll2009l.l201ll . herefater T09). The various progenitor models 
with varying mass, radii and explosion energy are summa¬ 
rized in Table [H and come fro m stellar evolutionary mod¬ 
els bv lUmeda fc Nomotol |200^. The first set of models are 
RSG with R/Rq = 500 — 1400 but we also add a second and 
a third set of models (not present originally in TOO) with a 
modification of the H/He envelope structure. The models 
in the second set have smaller radii of R/Rq = 80 — 400, 
and the models in the third set have dense envelopes with 



Figure 2. Spectral energy distributions at three characteristic 
times for a standard NSIO model of a red supergiant progenitor 
with 500/2©, 15M© and 1 foe. The mean wavelengths of filters 
NUVgalex and Sloan griz are shown. 

intermediate radii of R/Rq = 300 — 600 (see Figure [3|. All 
models assume a metallicity of Z = 0.02. Higher metallicity 
leads to more extended pre-SN structures, so its effect on 
the shock and rise of the light-curve is similar to progeni- 
tor radius. Explosive n ucleosynthesis follows the method of 
ITominaga et al.l 1 I 2 OOII) and radi ation hydrodynamics i s per- 
formed with the code STELLA (iBlinnikov et akll 19981 . l2000l . 
I 2 OO 6 I) . This model was successfully used in T09 for SNLS- 
04D2dc and its GALEX counterpart. 


3 DATA 

The data used in this analysis comes from two rolling 
searches designed primarily to discover type la supernovae 
at intermediate redshifts: the SDSS-SN with SN la redshifts 
of 0.1 < z < 0.4 and the SNLS with SN la redshifts of 
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Table 1. Progenitor Models from T09: main sequence mass, 
pre-explosion radius and explosion energy. A fixed metallicty of 
Z=0.02 is assumed. 


M/Mq 

R/Rq 

E/(10^^erg) 

15 

500 

1 

20 

800 

1 

25 

1200 

1 

25 

1200 

9 

25 

1200 

19 

30 

1400 

1 

13 

80 

0.9 

13 

80 

9.9 

13 

200 

0.9 

13 

200 

9.9 

13 

400 

0.9 

13 

400 

9.9 

13 

316 

1 

13 

398 

1 

13 

564 

1 



Figure 3. Density profiles of standard hydrodynamical models 
by T09 (pink), models with small radii (green) and dense models 
(black). 


0.1 < z < 1.1. These non-targeted searches have the ad¬ 
vantage of simultaneously discovering many SNe of differ¬ 
ent types and sampling light-curves at all epochs, including 
non-detections prior to explosion and the scarcely covered 
rise-times of SNe II. The non-la SNe represent an exquisite 
dataset which to-date has not been fully exp l ored (with a 


small number of excer 

)tions iTaddia et al. 20141: Tavlor et al. 

2014 

:lD’Andrea et al.ll2010l: 

Bazin et al. 20091: iNugent et al. 


2009) . 


1. SDSS-SN 


The SDSS-SN was one component of the SDSS-II projects 
using the 2.5m te l escope at Apache Point Observatory 
jGunn et al.l l200^ : I York et al.l[2000h . It repeatedly im- 
aged 300 deg^ every 3 -4 days in rest-frame in ugriz filters 
dFukugita et al.lll996l) w ith a spectroscopic fo l lowup pro¬ 
gram ( Sako et al. 2008h . In ID Andrea et H] ll201ol) . the 
light-curves of 34 spectroscopically confirmed SNe IIP in 
the redshift range 0.027 < z < 0.1 44 were published; 
and more recently, ISako etliD ll2014h presented the full 
dataset of more than 10000 transients discovered during 


the fall 2005-2007 campaigns. This full dataset comprises 
62 s pectroscopically confir med SNe II (including those 
from Id ’A ndrea et ahl l2010l ) but it also has many more 
photo metric SN II candidates with host galax y spectro¬ 
scopic dSako et al.1 20lJ: Olmst ead et al.ll2014l ) or photo¬ 
metric (DR8, lAihara et al.l 201ll) redshifts. 

2. SNLS __ 

The SNLS (IPerrett et al.ll201ol ) used the deep component 
of the 3.6m Canada-France-Hawaii Telescope (CFHT) 
Legacy survey to image four 1 deg^ fields repeatedly every 
2-3 days in rest-fr ame in filters QMrMiM ZM, similar to the 
SDSS passbands (iRegnault et al.ll2009l ). and with a spec- 
troscopic followup component at va r ious large telescopes 
( Howell_et_^ l2005l: [Hlis et al.l 20081: Bronder et ahllioosl : 


iBalland et al. l200^nWalker et al.l l201lj)^ A small frac¬ 

tion of the SNe IIP were used in Nugent et al.l (l2006l) 
for distance estimation analyses. The full dataset of pho¬ 
tometric transients consists of more than 6 000 objects, 
for which a fraction also have spectroscopi c (iDavis et al.l 
1200^ Le_^evreet_aljl200^ LilI^_etSTl200^ or photomet¬ 
ric ( Ilbert et al.ll200a : lFreeman et al.ll2009l ) redshifts from 
their host galaxies. 


We search for SN II candidates within the full pho¬ 
tometric sample of both surveys. For this, we use a sim¬ 
ple identification tool based on the shallowness of post- 
m aximum slopes in d ifferent bands (see Gonzalez method 
in iKessler et al.ll2010l) and a posterior visual inspection to 
eliminate false candidates. With this method we recover 
most spectroscopically classified SNe 11. However, given that 
the general spectroscopic classification of the SDSS-SN and 
SNLS may include other core-collapse SNe such as interact¬ 
ing SNe Iln and SNe Ilb, we revise all spectr oscopic classi¬ 
fications using SNID llBlondin fc Tonr'^l2007l ) and find that 
a 20-30% fraction of them (9 for the SDSS-SN and 3 for the 
SNLS) are not SNe IIP/L and were not classified so photo¬ 
metrically either. Nonetheless, we caution that there are a 
small number of events which are photometrically classified 
as SNe II, but spectroscopically are consistent with other 
SN types. This is a source of contamination and is further 
investigated in section [S] 

Our final SN II candidates are divided into a “golden 
sample” which need to fulfill following requirements: a) a 
spectroscopic redshift from the SN or the host, b) a spec¬ 
troscopic classification or a good photometric classification 
as a SN II and c) at least one data point between the last 
non-detection and the date of maximum in a given band. 
Such cuts result in 48 SNe for the SDSS-SN and 38 for the 
SNLS. A less restrictive “silver sample”, which includes the 
golden sample, requires a) a spectroscopic or photometric 
redshift, b) a spectroscopic classification or a photometric 
classification and c) at least one non-detection and an esti¬ 
mate of the date of maximum for a given band. As opposed 
to the golden sample, the silver sample allows photometric 
redshifts and does not necessarily require data within the 
rise. The silver sample contains 131 objects for the SDSS- 
SN in the redshift range 0.1 < z < 0.6 and 92 for the SNLS 
in the redshift range 0.1 < 2 < 0.8. A summary of these 
samples is presented in Table [2l The redshift distributions 
are shown in Figure |3] We also show a comparison of SN II 
candidates that have both, spectroscopic and photometric, 
redshifts. 
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Figure 4. Left: Redshift (spectroscopic and host photometric) distribution of silver samples for SDSS-SN (solid orange) and SNLS 
(diagonal blue) SNe 11. Golden samples are shown in dashed and dotted lines, respectively. Right: Comparison of photometric redshifts 
from the host with spectroscopic redshifts from the SN or the host for SDSS-SN (orange circles) and SNLS (blue squares) SN II candidates 
with spectroscopic information. Smaller open symbols are other transients from each survey. The solid line is a 1:1 relation and dotted 
lines represent 10% redshift error. 


Table 2. Number of SNe in the SDSS-SN and SNLS silver and 
golden samples. Spectroscopically confirmed SNe II in parenthe¬ 
sis. 


SDSS-SN 

SNLS 

silver 

golden 

silver 

golden 

131(24) 

48(24) 

92(14) 

38(14) 


in Aeff of the order of hundreds of A. Errors coming from the 
uncertainties on the observed colors are quite small, on the 
order of a few A. All these errors are added in quadrature 
and included in the analysis. 

4.2 Rise-times 


4 ANALYSIS TOOLS 

In this section, we present the different tools we use to 
analyze our photometric da ta. Instead of performing K- 
corrections llKim et al.lll996l) to correct light-curve data to 
the rest-frame, a process which has not been investigated 
enough for SNe II, we work directly in the observer-frame 
taking into accou nt the effective wave length. This is similar 
to the approach of iTaddia et al.l ll2014l) for stripped-envelope 
SNe. 

We work in flux (/) space and investigate each filter 
individually. Some example light-curves are shown in Fig¬ 
ure [5] For each SN and filter, we measure: a) the effective 
wavelength, b) the rise-time in each filter, and c) the linear 
slope (in magnitude space) after maximum. 


4.1 Effective wavelengths 

We measure effective wavelengths, Aefr, of each filter for 
each SN since the objects are found at varying redshift. 
We use the spectral energy distribution (SEP) template at 
maximum given by iNueent et al.l (l2002l n redshifted to the 
observer-frame and warped to m atch the observed colors at 
maximum (see lHsiao et al.ll2007l ). If we use observed spec¬ 
tra at e arly times instead (e.g . SN2007od 6 days aftex ex¬ 
plosion, iGutierrez et al.ll2015h . we obtain consistent effec¬ 
tive wavelengths to within 50A. Uncertainties in photomet¬ 
ric redshifts for some SNe in the silver sample result in errors 


To measure rise-times, we need to estimate a starting date 
and an end date of the rise for each band. For the starting 
date we identify two methods. The first consists of calculat¬ 
ing the mid-point between the epoch of last non-detection 
and the epoch of first detection. Although straightforward, 
this estimate can lead to systematic effects when obtaining 
non-detections with different limiting magnitudes. The sec¬ 
ond method improves the estimate on the explosion date 
by performing a simple power law fit to the data prior to 
maximum and including non-detections: 


f{t) = 


a{t - tLT)" 
0 


if t > 
if t < tL°p 


( 1 ) 


with the error obtained directly from the fit. The explosion 
date obtained in this way, tg°p, does not change significantly 
with the range of epochs used in the fit. The power index 
n, on the other hand, is affected by this limit and a lower 
upper range is preferred to properly estimate the early shape 
of the light-curve. To estimate the power-law, we therefore 
use data up to half of the flux at maximum. 

For the end date of the rise, we can use the date of 
maximum light in each band. Instead of directly taking the 
observed maximum which depends strongly on the cadence 
of the observed SN and generally leads to over-estimations 
in the case of short rise-times, we do a light- curve fit to 
the p henomenological function introduced by iBazin et al.l 
( 1200 ^ : 


^ https: //c3. Ibl. gov/nugent/nugent_teinplates.html 


fit) = A 


g —(i —io)/'^/aZZ 
\ ^o)/'t'rise 


+ B, 


( 2 ) 
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Figure 5. Observer-frame light-curve examples of SDSS-SN SN2005k b (left) and SNLS -06D3fr (right). Fluxes normalized to maximum 
of each filter are shown in different colors. Solid lines are fits to eq. [2] llBazin et al ] l200tl and dotted lines are fits of pre-maximum data 
to a power law. The small vertical ticks represent the fitted explosion dates, (lower ticks) and the end dates, (upper ticks) for 
each filter (some have been shifted by half a day for clarity). The arrows indicate the fitted maximum dates. The upper x-axis shows the 
restframe epochs with respect to the maxinmum in the bluest filter (u and g respectively). As an example, we quote the different rise-time 
measurements in g-band (corresponding to Aeg = 4511,3882A) for both SNe, SN 2005kb and SNLS-06D3fr: 4.4± 0.4,4.4± 1.2 and 

t;h“=4.7±0.4,4.5± 1.2 days; and in i-band (corresponding to Aeff= 7373,60291): t®pd^l0.1±1.7, 6.9±1.2 and tjhjf=41.4±1.3, 7.9±1.2. 


where A, B, to, Tfaii and r^ise are free parameters. We can 
get a date of maximum, tmax) from the derivative of this 
function: 

Cax =t0+ Tris. log ( , (3) 

\Trise + Tfall / 

and an error propagated from the la uncertainties on the 
parameters from the fit. 

In some cases, the date of maximum flux in the red¬ 
der bands can occur during the plateau, much later than 
the end of the rise and much later than the maximum in 
bluer bands (see Figure [3. Although the time between ex¬ 
plosion and maximum light is easily estimated, and has been 
used in the past to measure rise-times of nearby SNe II, it 
can be strongly affected by other physical processes like hy¬ 
drogen recombination at longer wavelengths. To avoid over¬ 
estimating rise-times, we use the derivative of the same phe¬ 
nomenological function and find the epoch at which it be- 
comes lower than O.Olmag/day, in a similar fashion to 
I Gall et ^ ll2015l ). This end date is similar to the date of 
maximum at shorter wavelengths but avoids possible rising 
plateaus in redder bands. The error is calculated by prop¬ 
agating the errors on the fit parameters. Some examples of 
these estimates and fits are shown in Figure [S] We thus 
present rise-time measurements based on the fitted explo¬ 
sion date, , and the end date from the Bazin function, 

> Baz, 

^end • 

■ end _ ,Baz ,pow r a\ 

^rise ^end ^exp • 

For comparison with previous works, we also show the 
rise-times based on the maximum date: 

.max _ ,Baz ,pow /r\ 

^rise ^max ^exp • 

As will be shown, our main conclusions are not affected 
by the use of different measurements of the starting date 
end dates for the rise-time calculation. 


4.3 Post-maximum slopes 

We measure a post-maximum slope in each filter by doing 
a linear fit in magnitudes to the data after maximum up to 
80 days past maximum. We require at least three epochs to 
ensure a proper ht. This gives a general quantitative idea of 
the shape of the light-curv e after maximum . It is r eminiscent 
of the “s2” plateau slope of lAnderson et akl (12014 ) except we 
do not fit a primary “si” decline. 


5 RESULTS 

5.1 Rise-times versus effective wavelength 

We present the rise-times as a function of effective wave¬ 
length for the SDSS-SN and SNLS, and compare them to 
analytical and hydrodynamical models. For consistency, we 
calculate the rise-times for data and models in the same 
manner. Figures |^and [7]show the rise-times calculated with 
and t®ax) respectively, as a function of effective wave¬ 
length. We show in colors the observer-frame filters from 
which we measured each rise-time, so that typically a SN 
with good coverage has 4 (SNLS) or 5 (SDSS-SN) values in 
these plots, each at different wavelength. The observed color 
of the SN influences the effective wavelength but it is mostly 
the redshift that gives rise to the diversity for a given photo¬ 
metric band. Since the SNLS spans a larger redshift range, 
the dispersion on Aefi for a given filter is larger and one can 
see overlap between two filters for the same wavelength. A 
summary of the rise-times is presented in Table [21 

The bulk of the rise-times, trSe, lies below 10 days 
for wavelengths less than 4000A and increases to about 
12-15 days beyond 7000A. If one uses the maximum rise- 
times, they are of the order of 15-20 days be¬ 

yond 7000A. The population of SNe II clearly has steep 
rise-times, which are signihcantly faster than those es- 
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SDSS-SN 


SNLS 



Effective wavelength [A] 



Figure 6. Calculated rise-times, for the SDSS-SN (left) and SNLS (right) samples as a function of effective wavelength. Each point 
is the rise-time of a SN at a given band. All SNe from the silver sample are shown and the SNe from the golden sample have black 
contours. Different colors indicate the observed filter and larger symbols denote the median values for characteristic wavelengths. Some 
representative hydrodynamical models by T09 are shown: standard models of large radii in pink, models of small radii in green and dense 
models in black. 


Table 3. Median and deviation on the median for two rise-times measurements, and for the SDSS-SN and SNLS silver and 

golden samples at different effective wavelengths. 


Effective wavelength 


SDSS-SN 



SNLS 


(A) 

silver 

golden 

silver 

golden 


e^(d) 

i™(d) 

i“^(d) 
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Figure 7. Calculated maximum rise-times, for the samples of SDSS-SN (left) and SNLS (right) as a function of effective wavelength. 
Each point is the rise-time of a SN at a given photometric band. Black contours indicate SNe from the golden samples, the rest are from 
the silver samples. Different colors indicate the observed filter and larger symbols denote the median values for characteristic wavelengths. 
Some representative analytical models by NSIO are shown for comparison: a RSG progenitor with ISMq, 500Rq and Ifoe explosion 
energy (solid black line), varying the mass to 10 or 25Mq (dashed and dot-dashed maroon), or varying the radius to 50 or 2500Rq 
(dashed and dot-dashed blue), or varying the energy to 0.1 or 10 foe (dashed and dot-dashed green). 
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timated for stripp ed-envelope SNe (12-23 days in r/R 


iTaddia et alJ 

20141) or standard thermonuclear SNe la (16- 

25 days in B, 

Firth et al.ll2ni4l:lGonzalez-Gaitan et al.ll2012l: 

iHavden et al. 

l2010h. This is a possible indication of pure 


shock cooling of an extended progenitor star powering the 
early light-curve of SNe II, as opposed to compact stripped- 
envelope SNe or SNe la which are generally caught when 
the light-curve is already powered by radioactivity, e.g. the 
shock cooling for W olf-Rayet stars is aronnd 100 seconds 


dKistler et al. l2013l l as opposed to the characteristic de¬ 
cay time of 8.8 days for ®®Ni. Previous works of nearby 
SNe have already shown that individual SNe II have sharp 
rises t o maximum: SN 1999 em rises in 6, 8 and 10 days in 
UBV \heonaxA et al.l 200 2^. SN 1999gi rises in 8 and 12 


days in i?I/| Leonard et ^l2002al~). SN 2004et r ises in 9, 10 


16, 21 and 25 days in UBVRI i Sahu et'^l2006h . SN 2012aw 


rises in less than 3 days in UV and in 8, 11, 15, 22 and 24 
days in UBVRI and SN2012A rises in 10, 12, 17 and 24 
days in BVRI llXomasella et al.ll2013ll . All these values are 
consistent with our measurements. 

In the figures we also show several analytical models 
from NSIO and hydrodynamical models from TOO. We see 
that for both, SDSS-SN and SNLS, the standard models of 
TOO from RSG progenitors with radii 5OO-14OO7?0 predict 
long rise-times and therefore lie consistently above the mea¬ 
sured rise-times (pink lines in Figure O. In the case of a 
second set of new hydrodynamical models with small radii, 
80-400 7?0 (green lines), we see that the predicted rise-times 
are in better agreement with the observed rise-times. This 
is clear evidence for smaller radii than typical assumed RSG 
radii. 

With the analytical models, we see a similar prediction: 
the rise-times of a 500 Rq progenitor are roughly consistent 
with the SDSS-SN and SNLS maximum rise-times for NSIO 
(see Figure [7|). Furthermore, we show other models vary¬ 
ing the mass, radius and explosion energy and find that the 
largest variation in rise-time duration can only be accounted 
for by radius. As a matter of fact, the luminosity and tem¬ 
perature in their model equations scale more strongly with 
radius than with other parameters. 


5.2 Distribution of rise-times 

We find that there is a distribution of rise-times for a given 
wavelength that is larger than the measured rise-time uncer¬ 
tainties. At Aefi<4000A, the deviation on the median (MAD) 
is about 4 days. This indicat es that besides their b rightness 
and post-maximum decline (I Anderson et al.ll2014fl . SNe II 
constitute a heterogeneous sample during the rising part of 
their light-curve as well. To further study this, we interpo¬ 
late the rise-times (trise) of each SN to the common effective 
wavelength of the SDSS p-band at rest-frame, i.e. 4722A, 
and plot the distribution in FignreO One can see that the 
majority of SNe II have rise-times below 10 days at this 
wavelength although a population of objects with long rise- 
times exists and is investigated in next section. To test how 
our large uncertainties in the rise-times affect this distribn- 
tion, we perform a Monte Garlo analysis of 100 realizations 
in which each rise-time is randomly shifted according to its 
one (7 uncertainty. We then calculate the mean of all simula¬ 
tions to obtain a characteristic rise-time of 7.5±0.4 (6.0±0.3) 
and 10.3 ± 0.4 (10.7 ± 0.7) days for the silver (golden) sam- 



SNLS 

SDSS-SN ■ 


Figure 8. Histogram of rise-times at rest-frame g-band 

(Aefj=4722A) for the silver samples of the SDSS-SN (solid or¬ 
ange) and SNLS (diagonal blue). The golden samples are shown in 
dashed and dotted lines, respectively. Median MC values are indi¬ 
cated with ticks at the top: 7.5'Gie and 6 . 0^2 * days for the silver 
and golden samples of the SDSS, and 10 . 3 ^ 4 '| and ® days 

for the silver and golden samples of the SNLS, where the quoted 
uncertainties contain 84% of the distribution. The vertical dotted 
line indicates the limit of long rise-time SNe (see section ISAjl. The 
horizontal error bar represent the typical median error on the rise 
of 2.5 (SNLS) and 1.6 (SDSS-SN) days. 


pies of the SDSS-SN and SNLS, respectively. If we take both 
surveys together we obtain 8.9 ± 0.2 and 7.5 ± 0.3 days for 
silver and golden samples. We also calculate for each simu¬ 
lation the lower and upper rise-times for which 84% of the 
population is contained and obtain a final mean range of 4.5- 
16.6 days. As mentioned in the previous section, all nearby 
SNe have rise-times of 8-11 days in R-band (Aeff= 439lA), 
in good agreement with this distribution. 

The large discrepancy between both surveys is found in¬ 
dependently of the way we calculate the rise-time. It is due 
to the difference in cadence for both surveys. In hgure|^we 
show the measured cadence during the rise for each SN light- 
curve in our samples. The cadence is the mean rest-frame 
duration between data points calculated during the rise for 
a given SN light-curve. We obtain very different average ca¬ 
dences of 4.1 and 8.2 days for the SDSS-SN and SNLS. The 
mean cadence of the SNLS is of the order of the rise-time 
measure itself. Since we require at least one point during 
the rise to calculate the time of explosion in the power-law 
eq. m there will be a lack of objects with short rise-times, as 
is clearly seen in Figure |8] If we use the midpoint between 
last non-detection and the hrst detection to estimate the ex¬ 
plosion date instead, we are able to use more SNe, but the 
rise-times for the SNLS are still consistenly larger due to 
the long time range without any detection (see simulation 
in appendix IQ. Thus, we believe that the rise-time distri¬ 
bution of the SDSS-SN is more complete and we stress the 
possibility that there might be even shorter rise-times that 
are not probed with the cadence of the SDSS-SN. The main 
conclusions of possible small radiie are not affected by this 
systematic; in fact they are strengthened. 
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Figure 9. Cadence as a function of effective wavelength for the SDSS-SN (left) and SNLS (right). Each point is the rise cadence a SN 
at a given filter. SNe from the golden sample have black contours. Different colors indicate the observed filter and larger symbols denote 
the median values for characteristic wavelengths. 


5.3 Long rise-time SNe 

Despite the fast rise-times found for most of our SNe, we find 
a fraction of SNe that have much longer rise-times, i.e. larger 
than 15 days (see Figure ITOll . Such different rise durations 
could indicate that these SNe are not really normal SNe II 
but are powered by different energy sour ces like radioactiv¬ 
ity for peculiar SNe II such as SN 1987A l|Utrobin fc Chug^ 
1201 ll : IPastorello et al.ll2ni2h or interaction for SNe Iln, or 
even mo re exotic events suc h as the super-luminous SNLS- 
07D2bv l|Howell et aLlfioI^ . Some of these different types 
may have shallow slopes like normal SNe II and could be se¬ 
lected with our classification technique that is based only on 
post-maximum photometry. To investigate this further, we 
take a closer look at the SNe originally classified as SNe II 
based on post-maximum slopes that have y-band rest-frame 
rise-times at least one sigma above 20 days. We indeed find 
that one of these was rejected based on its spectrum being 
consistent with SN Iln showing narrow absorption Ha lines 
(SN2006ix). Even though it was classified as SN II photo¬ 
metrically, its long rise-time actually reveals that it is not 
a standard SN II. Of all spectroscopically confirmed long 
rise-time SNe, the only one that is clearly a SN II is SNLS- 
07D2an, but the spectrum, particularly the Ha P-Cygni pro¬ 
file, is much more consistent with that of peculiar SN II 
SN 1987A. We note that this SN has much shallower post¬ 
maximum slopes, characteristic of SNe IIP, as opposed to 
other SN 1987A-like which decay faster. 

This provides strong evidence that normal SNe II are 
constrained to short rise-times and other possible hydrogen- 
rich SNe have longer rise-times. These SNe were identified 
based on their exceptionally long rise-times besides their 
spectra, demonstrating the power of photometric classifica¬ 
tion from just the early part of the light-curve. If we se¬ 
lect out all photometrically identified SNe II that have long 
rise-times of a sigma more than 15 days, we find that a 
7-10% (7% for SDSS-SN and 10% for SNLS) has a longer 
pre-maximum behaviour from normal SNe II. We argue that 
most of these SNe are actually of a different origin but we 
keep them in the analysis since we do not have a spectro¬ 


scopic confirmation. If we were to take them out when calcu¬ 
lating the median values shown in the figures and in Table |31 
we would obtain median Monte Carlo rise-times shorter by 
day. 


5.4 Rise-time differences versus effective 
wavelength 

One can see in Figures|n]and [7]that the rise-times are shorter 
at bluer wavelengths, i.e. bluer bands peak before redder 
bands, which is expected from cooling of the photospheric 
temperature shifting the peak of the emission towards longer 
wavelengths as time evolves. The relation can be further seen 
in the rise-time differences of each SN respect to an interpo¬ 
lated rise-time at rest-frame g-band (Figure [TT}. This trend 
is similar to the behavi our observed in oth er SNe such as 
stripped-envelope SNe (iTaddia et al.l [20l3l . We fit a sec¬ 
ond order polynomial aX^ -|- &A -|- c to these differences and 
find that they represent the data better than a linear trend 
according to an F-test. Comparing these fits with stripped- 
envelope SNe, we do not find a statistically significant differ¬ 
ence although the physical processes driving the rise are of 
a different nature: for stripped-envelope SNe the rise-times 
are powered by radioactive nickel heating whereas SNe II 
mostly follow the shock-heated cooling at these epochs. 

Although the models also show an increase in rise-time 
at longer wavelength, the relation is generally steeper for 
the analytical models than for the observations, except for 
very small radii (~ 50 R©), typical of BSC progenitors, that 
are inconsistent with the bulk of the absolute measured rise- 
times. The set of hydrodynamical models of progenitors with 
small radii 2007?©), on the other hand, are more con¬ 
sistent with the differences between rise-times at different 
wavelengths. Regardless, both analytical and hydrodynam¬ 
ical models favor again small radii based on the rise-time 
differences at varying wavelengths. We have compared the 
photospheric temperature evolution of models with differ¬ 
ent radii and find that larger pre-SN radii lead to a slower 
temperature evolution. The slower temperature evolution re- 
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Figure 10. Observer-frame light-curve examples of two long rise-time SNe: SDSS-SN SN2006ix, possibly a SN Iln (left), and SNLS 
07D2an, a SN 1 987A-like SN II (ri ght). Fluxes normalized to maximum of each filter are shown in different colors. Solid lines are fits 
to equation 1 of iBazin et al.l ll2009lV The upper 3::-axis shows rest-frame epochs with respect to maximum in the bluest band {u and g 
respectively). 


Table 4. Median and deviation on the median for the difference between the rise-times, different wavelengths respect to the 

effective wavelength at rest-frame p-band (4722A) for the SDSS-SN and SNLS silver and golden samples. 
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Figure 11. Difference in rise-times, respect to an interpolated rise-time at ^-band rest-frame (Aeff=4722A) for the SDSS-SN (left) 

and SNLS (right) samples. Points with black contours are the golden samples and larger symbols are the median values at characteristic 
wavelengths. The best second-order polynomial fit is shown as a solid line (orange for SDSS-SN and blue for SNLS). Overplotted are 
representative analytical models from RWIO for extreme radii of SORq (blue dotted) and 2500770 (blue dashed), and hydrodynamical 
models from T09 for a standard SOORq (pink), a more compact 200770 (green) and a dense progenitor with 316770 (black). 
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Figure 12. Aggregate light-curves of 11 SNe II from the SDSS- 
SN and SNLS that have fitted power-law rise-times between 6.5 
and 7.5 days. The solid lines shows power-laws of n = 1.5 (thick) 
and n = 1.0 (thin). The dotted line is the hydrodynamical model 
by T09 with R=200 Rq (with n = 1.3) while the dashed and 
dot-dashed are the analytical models by NSIO and RWll with 
R=400 Rq (with n = 1.0 and n = 0.8 respectively). 

suits in later peaks at red wavelengths, leading to a steeper 
dependence of rise-time with wavelength. 

5.5 Power-law 

The power-law of the rising light-curves can provide im- 
portant informatio n on the explosion physics of SNe (e.g. 
IPiro fc Nakarll201,3l ). From the fits to equation [T] we can ob¬ 
tain the power law index n for SNe with sufficient early 
coverage. The result of the fit is highly dependent on the 
range of epochs used in the fit. For simplicity, we take data 
before half the maximum flux in a given band. As the rise¬ 
time of SNe II is fast, this is only possible for a small subset 
of SNe II. We obtain a median power index of O.QGlJ y® for 
the SDSS-SN and 0.9ljlJ'5g for the SNLS, where the up¬ 
per and lower uncertainties represent the one a dispersion 
of the distribution. Although most of the SNe that we can 
calculate this for require data in the rise and are therefore 
biased to longer rise-times, we do not find any evidence for 
a trend of the power-law with effective wavelength or rise¬ 
time duration, and this is confirmed for the analytical and 
hydrodynamical models. Taking advantage of the combined 
statistics of multiple light-curves, we can make aggregate 
light-curves for SNe of a given rise-time duration to calcu¬ 
late the power-law. In figure fT^ we add the light-curves of 
11 SNe with rise-times between 6.5 and 7.5 days to a com¬ 
mon date of explosion obtained through individual power- 
law fits. We can see some scatter in the rise which may come 
from uncertainties in the explosion dates of each SN. Nev¬ 
ertheless, we obtain consistently low power-law indexes of 
< 1.4 for several aggregate light-curves of different rise-time 
durations between 5 and 15 days. 

We note that similar power-law fits to both, analyti¬ 
cal and hydrodynamical models, also predict such shallow 
power-laws. This is n ot surprising: based on the analyti¬ 
cal equations of NSIO, IPiro fc Nakarj (l2013l ) show that the 
UV and optical luminosity of a light-curve dominated by 


shock-heated cooling rises as T'^. When a light-curve rises 
more steeply than this, it probably means that other effects 
such as radioactive heating or recombination have started 
to play a role in the energy input. Given that most of our 
individual SNe and all aggregate light-curves have power- 
law indexes lower than 1.5, this provides strong evidence 
that the early light-curves of most SNe II are powered by 
shock-heated cooling. For ah models, we see that the post¬ 
maximum slopes decline much faster than for the observa¬ 
tions. This will be further discussed in section [5Tl 

5.6 Inferred pre-explosion radii 

If the dispersion in the rise-times found in previous sections 
is real, the most plausible progenitor property in the analyti¬ 
cal models that can account for this is radius (see Figure 0. 
Pursuing this assumption further, one can then use NSIO 
and RWll to infer radii based solely on the rise-time dura¬ 
tion and keeping the other parameters fixed at 15 Mq and 
1 foe. Taking the multiple fitted rise-times at different wave¬ 
lengths of each SN weighted by their error, we can then com¬ 
pare to the values measured for analytical models and ob¬ 
tain an optimal radius. The cumulative distribution of these 
radii for the golden samples of both SDSS-SN and SNLS 
together are shown in Figure [13] for the two models. The 
distributions predicted by NSIO and RWll are consistent 
(Kolmogorov-Smirnov KS-test probability of being drawn 
from the same population of 0.53) and the median radii are 
remarkably similar: 3231 jgg and 336ljgg7?0 (SShlgio 
4181253^0 silver samples). The uncertainties here 

quoted represent the 84% dispersion range. If we perform a 
Monte Carlo simulation in which we measure new radii from 
simulated rise-times based on the uncertainties and covari¬ 
ances, we find consistent mean radii of all simulations of 
3381202 and 34712®° Re (39512^ and 4251 26^7?© for the 
silver samples) for each analytical model respectively. 

We find that less than 5% of the objects are consis¬ 
tent with BSG progenitors powered exclusively by shock 
cooling of radii less than 1007?©, although we are limited 
by the cadence of the surveys. This is roughly in agree- 
ment with the 2% fraction inferred for SN 1987A events 
llKleiser et akllioill ). On the other hand, some objects hav¬ 
ing long rise-times have radii of several thousands Rq. As 
explained in section ESI these objects probably have a dif¬ 
ferent powering source. If we do not take into account the 
identified peculiar SNe with long rise-times, we also find less 
than 10% of our sample having longer radii than 12007?©. 
These radii are also confirmed qualitatively with the hy¬ 
drodynamical models of TOO (see appendix El for non-LTE 
treatment). In Figure[6]one can see that their models of radii 
much larger than 5007?© are ruled out and in fact smaller 
radii of ~ 100 — 4007?© are favoured. Small radii have been 
recently inferred through modelling o f individual sub lumi- 
nous SNe II like S N 2008in (1267?©, IRov et al.l[201ll) and 
SN 2009N (2877?r7.. lTakats et al.l[2^) . 

Typical RSG span radii of 1 00-16007?© llLevesaue et al.l 
I2OO5I : lArrovo- Torres et al.l l2013l ) and interferometric stud¬ 
ies confirm the existence of some very extended nearby 
stars with more th an 9007?© llHaubois et al.l I 2009 I : 
IWittkowski et ah! |2012| ). According to our findings, most 
SNe II are produced by RSGs with radii at the lower end 
of this distribution. This is shown in Figure [13] where we 
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Figure 13. Cumulative distribution of inferred radii for the 
golden samples of SDSS-SN and SNLS together calculated from 
multi-wavelength rise-times, with models NSIO (solid gray) 

and RWll (diagonal red). Vertical dashed lines show the median 
radii for both models. The dashe d cyan line shows the distribu¬ 
tion of RSG radii from the MW llLevesgue et al.ll2005h and the 
dashed green line the distribution of RSG radii from the MC 
jLevesQue et al.ll2006h . Dotted lines are Salpeter IMF distribu¬ 
tions of expected radii assuming R ~ with initial masses 

of 8-30Mq. Median errors on the radii are shown as horizontal 
bars in the right. 


present the distribution of RSG radi i for the Milky Way 
(MW ) and Magellanic Clouds (MC) (|Levesaue et alJ[20Qa . 
l2006fl . These were obtained with Stefan Boltzman’s law 
from the effective temperatures and bolometric luminosi- 
ties inferred with MARCS stellar atm osphere models (see 
iPlez et af]ll992l : rGustafsson et al.lll975l l applied on observed 
spectrophotometric data. The distributions are inconsistent 
with the radii inferred from our SN rise-time analysis (KS- 
test of being drawn from the same population of less than 
10“® and 10“^® for the MW and MC respectively). Since 
smaller stars are also less massive, one is tempted to as¬ 
sociate this observation with the seemingly lack of massive 
RSG progenitor of more than ~ 18 Mq in pre-explosion im¬ 
ages JSmartt et al.l[2009l ~). Figure [14] shows RSG radii against 
the masses obtain ed from the mass-luminosity relation for 
the MW and MC dLevesaue et ahilioosl . I2OO6II . We include 
our one u (84%) upper bounds on radii for NSIO and RWll 
models. They roughly correspond to upper limits of 19 Mq, 
consistent with the limits from pre-explosions masses. This 
would be an independent confirmation of their result. 

Nevertheless, the disagreement between our inferred 
radii and the observed RSG from the MW could just be 
a result of the incompleteness of their survey for which the 
smaller RSG are not all observed and are therefore missing 
in the distributions. Instead, we can do a very basic com¬ 
parison with the expectations of a standard Salpeter Initial 
Mass Function (IMF) assuming all SNe II arise from massive 
single stars between 8 and 30Mq (the upper limit has very 
little impact here). To obtain radii, we use a mass-radius re¬ 
lation of R/Rq = 1A{M/Mq)^'^ from a log-log linear fit to 
the MW RSG of iLevesaue et al.l (l2006tl . The resulting sim¬ 
plistic radius distributions of this IMF are shown in dotted 
lines in Figure [T^ The normalization is obtained by fitting 
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Figure 14. Radius versus mass for RSG in the MW and MC 
llLevesgue et al.ll2005l . I 2 OO 6 I ) and 84% upper limits of our radius 
distribution with NSIO and RWll models as solid lines. Dashed 
lines are the approximate mass correspondance for the MW. 

the IMF to the NSIO and RWll distributions. In this case, 
our inferred distributions are more consistent with the sim¬ 
ple IMF distributions, especially for RWll (KS-test of being 
drawn from the same population of 0.03-0.09 for RWll and 
0.09-0.2 for NSIO). 

On the other hand, another possible reconciliation 
stems from recent RSG observations extending to the IR 
which suggest that they are hotter than previously thought 
llDavies et al.l[201^ . which can re sult in a reduction b y up 
to 30% of the inferred RSG radii (IPessart et al.l[2^013ll . Ar¬ 
bitrarily shifting the distributions of RSG radii in the MW 
by this amount results in a population that is more consis¬ 
tent with our inferred radii (KS-test of being drawn from the 
same population of 0.02 — 0.1). Furthermore, there are many 
uncertainties regarding the stellar models used to obtain 
these radii, in particular the mixing length can greatly affect 
the i nferred radius (e.g. IPeng fc Xiona[200ll : iMevnet et al.l 
l2014h . These difficulties raise the question if RSG radii are 
systematically over-estimated and the distribution is more 
similar to what we find based on SN rise-times. 

Finally, metallicity can affect the rise-t imes through 
changes in the pre-SN progenitor struc ture (iLangerl Il99ll : 
iGhieffi et al.ll20(^ : iTominaga et al.ll201ll l. Low er metallicity 
result s in shorter radii: according to models bv lGhiefh et al.l 
1 I 2 OO 3 II a reduction from solar metallicity Z = 0.02 to 
Z — 0.001 reduces the pre-SN radius by Ai? = 150 — GOORq 
for initial masses of 13-15M0. This could explain our find¬ 
ings but it would mean that most SN II are hosted in low 
metallicity environments which i s not a general observ ation 
for nearby targeted hosts (e.g. lAnderson et al.l l2010l i but 
should be tested for our sample. The distribution of RSGs 
in the low-metallicity MG does not support this either, nor 
the longer rise-times found for the SNLS, which at higher 
redshift is expected to have lower mean metallicity. 

5.7 Slopes versus effective wavelength 

Here we analyze the post-maximum linear slopes measured 
in each band as a function of effective wavelength for both 
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SDSS-SN SNLS 




Effective wavelength [A] Effective wavelength [A] 

Figure 15. Post-maximum slopes as a function of effective wavelength for the SDSS-SN (left) and SNLS (right) silver (all points) and 
golden (points with black contours). Standard (pink), small-radii (green) and dense (black) hydrodynamical models from T09 are shown 
in lines. We do not present analytical models because NSIO and the simple model of RWll that we use do not treat the recombination 
of hydrogen. 


samples (see figure nsi). At the shortest wavelengths, the 
light-curves decay much faster than at longer wavelengths, a 
consequence of the temperature evolution of SNe II. The de¬ 
cay rate of optical wavelengths above approximately SOOOA 
are more similar to each other reflecting the known hydro¬ 
gen recombination wave as it recedes through the ejecta. 
The slope and duration of this plateau phase are possibly 
affected by the mass of hydrogen present in the ejecta. This 
can be best seen in the hydrodynamical models of T09: 
models with large radii (pink) predict the right range of 
slopes whereas the set of models with small radii (green) 
overpredict the slopes. This result can be explained by 
a different ejecta mass or radius in the RSG progenitor 
jLitvinova fc NadezhinlllQSSl : |Popovlll993l l. RSG with large 
radii have a larger hydrogen mass that produces a shallower 
(and possibly longer) plateau and RSG with small radii have 
much less hydrogen and hence the recombination wave re¬ 
cedes too fast and does not produce a plateau. Co mparing 
our slopes with the results of lAnderson et ahl (l2014l l. we ob¬ 
tain a very consistent range of 0.3 — 4 magnitudes per 100 
days in the wavelength range of 5000-6000A (roughly equiv¬ 
alent to V-band), indicating that the populations at high-z 
are quite similar to \aw-z. We note that our slopes are a 
combination of their si and S 2 values. 

Together with the results presented in the previous sec¬ 
tions, we then require RSG progenitors of small radii to pro¬ 
duce short rise-times but with large hydrogen mass to pro¬ 
duce shallow post-maximum slopes. Thus a dense hydrogen 
envelope is needed which can be achieved through stronger 
convection in RSG envelopes that leads to more mixing. 
Indeed, there is discussion on t he variation of con vective 
mixing length in different stars dFerraro et_ah 2006ll which 
is po orly constrained for massive stars i Meakin &: ArnettI 
I2OO7I). Varying the mixing length severely affects the ra¬ 
dius dMaeder fc Mevn^ll987l ') : iDessart et al.l d2013h showed 
that varying i t from a = 1.1 to 3 in MESA STAR models 
dPaxton et al.ll201ll l changes the radius from 1100 to 500 
7 ? 0 , and it does not affect the He core characteristics since 
the envelopes are only influenced at the end of their lives. 


They also find that a small radius is needed to solve SN color 
evolution discrepancies of SN 1999em from stellar evolution¬ 
ary models coupled with SN radiative transfer. This opens 
up an interesting test that can be performed in the future, 
both on the current and other samples, by looking at the 
color evolution and its relation to rise-time. Also, different 
treatments of the mixing length theory for certain RSG en¬ 
velopes in which convective speed bec omes supersonic, pre¬ 
dict smaller radii for the same mass dPeng fc Xiond [2001! : 
Ijosselin fc Plejl200lfl . We note that enhanced rotation or 
core overshooting may also alter the ra dius but in the op¬ 
posite direction to more extended RSG dMaeder fc Mevnea 
I2OOOI: iDessart et al.l[20ill . 

We have included three additional hydrodynamical 
models with dense envelopes to test this hypothesis further. 
These models are shown as black lines in previous Figures. 
In Figure 1161 we directly compare the pre-maximum and 
post-maximum behaviors, i.e. rise-time and plateau decline, 
for both samples in light of different hydrodynamical mod¬ 
els. We do not see any bimodality indicative of two groups, 
SNe IIP and SNe IIL. There is however a trend for slower 
decliners to have longer rise-times and faster decliners to rise 
quicker (Spearman correlation of -0.39 and -0.27 for golden 
and silver samples). Dividing the population by the median 
rise-time (8.9 days), the distributions of slopes of SNe II 
with short rise-times is significantly different from the one 
with longer rise-times (KS-test of 0.002 for golden and sil¬ 
ver samples). In the traditional nomenclature, this means 
that SNe IIL tend to have steeper p ost-maximum slopes and 
shorter rise-times. iGall et al.l (l2015fl find a similar trend. The 
most straightforward interpretation based on the previous 
discussion is that SNe II that decline more rapidly during 
the plateau have smaller radii. If a smaller radius means less 
hydrogen, then this is naturally explained. 

The hydrodynamical models indeed predict a relation 
between the rise-time duration and the post-maximum de¬ 
cline. However, none of them is able to reproduce the 
parameter space in Fi gure JT6L This leads us to suggest 
an alternative scenario. iMoriva et al.l 1 I 2 OIII . Figure 7) and 
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I Baklanov et alJ (l2013l . Figure 1) present models for a family 
of SNe Iln, i.e. SNe IIP with circumstellar material (CSM) 
interaction. From their figures one can clearly see that the 
optical rise-times of their models are much shorter than pure 
SN IIP models; in fact they can often have values below 10 
days as we hud in this study. In those models the shock 
breakout occurs at larger distances, beyond the stellar sur¬ 
face, and the signal is therefore smeared in the CSM and its 
duration is prolonged. Consequently, the rise-time directly 
corresponds to shock breakout and not shock cooling. Would 
it then be possible that normal SNe II indeed have pre-SN 
mass-loss high enough to affect the shock and the rise-time 
but not so strong as to create narrow emissi on lines in the 
post-maximum spectra? ICezari et al.l ll2015ll recently con- 
hrmed a much quicker rise in NUV for SN IIP PSl-13arp 
and propose a similar explanation. If this were the case, 
then faster decliners (SNe IIL), that seem to have shorter 
rise-times, may have a larger or denser CSM interaction. 
The interaction would naturally explain their higher bright¬ 
ness as the kinetic energy is more efficiently con verted to 
thermal energ y in t he CSM (iMoriva et al.l uOlll l. In fact, 
IValenti et al.l 1 20151 ') recently detect moderate CSM inter¬ 
action for type IIL SN 2013by and suggest such a com- 
mon feature for al l SNe IIL. Moreover, the recent work by 
ISmith et al.l (l2015l l shows the existence of transitional ob¬ 
jects between classical SNe II and interacting SNe Iln such as 
PTFlliqb. These objects initially behave like SNe Iln with 
narrow emission lines that later disappear resembling more 
SNe 11. This is further indication that SNe II and SNe Iln 
may all be part of a continuous family with different CSM 
quantity or density in following increasing order: IIP-IIL-IIn. 
We note that in principle any SN type can be hidden behind 
a thick CSM, making it look like a SN Iln. lOfek et al.l (l2014l l 
actually find a large diversity of rise-times for SNe Iln. Based 
on our hndings, it is tempting to argue that only SNe Iln 
with short rise-times are SNe II inside a thick CSM, where 
the energy input is more strongly attributed to the shock 
breakout. 

Finally, a differ ent observation of interest by 
lAnderson et al.l ll20l4 l in which SN II P-band light- 
curves (but also observable in others bands such as B) are 
found to be described best with two declining slopes after 
maximum, si and S2, suggests that the first steeper decliner 
stems from exten ded cooling, which is related to the radius 
of the progenitor (iBersten et al.ll2012h . Such late transition 
times between si and S 2 require very large radii which, 
together with the fast rise-times we measure, can only be 
reconciled in this CSM scenario. 

All these observations point towards an important 
paradigm shift in which SNe II (IIP and IIL) also have 
pre-SN wind mass-loss of 10“^ — lower than 

SNe Iln of 0.01 — 0.1MQj/r“^. The obvious question arises 
on the coincidence of mass-loss and SN explosion. Ei¬ 
ther episodic mass-loss is extremely com mon or there is 
some m echanism that makes them coincide . iMnley &: Sokeil 
iI 20141 : ISokerl (120131 '): iQuataert fc Shiod3 1 20121 ') present 
such a possibility: the RSG envelope expands shortly be¬ 
fore explosion due to powerful nuclear burning and so 
called p-waves carrying nuclear energy from the core to 
the outer layers thro ugh core convection. This can lead to 
extended envelopes jMclev k, SokeJ I20l4 ISokeil l2013l j or 
mass-loss jQuataert fc Shiodd 20121 ') . which in both cases 



Figure 16. Post-maximum slopes versus rise-time for the silver 
samples of SDSS-SN (orange circles) and SNLS (blue squares). 
Golden samples have black contours. Measurements in all photo¬ 
metric bands are included. Selected standard (pink), small-radii 
(green) and dense (black) hydrodynamical models from T09 are 
shown in lines. 


results in shorter rise-times. Binary interaction could play 
an important role in recurrent mass-loss episodes and 
asymmetric CSM, as observed for some in teracting SNe 
jMauerhan fc Smithll20l3 : ISmith et al.ll2015l ). 

Instead of CSM, a very extended envelope can also 
produce short rise durations. As a matter of fact, ex- 
tended molecular atmospheres have been obse r ved fo r RSG 
dWittkowski et all l2012f : lArrovo- Torres et al.l l2013l . 120151) 
and densities in t he outer layers of RSG are larger than con¬ 
vection predicts dOhnaka et al.l[201ll . l2013l ). These explana¬ 
tions naturally fit the short rise-times, shallow plateaus and 
late transition of si to S2. 


6 SUMMARY 


We have measured the rise-times as a function of effective 
wavelength for a large sample of SNe II spanning a wide 
redshift range of 0 < 2 < 0.7 from the rolling surveys 
of SDSS-SN and SNLS. The rise-times increase with wave¬ 
length as expected for shock-heated cooling. The power-law 
indexes of the rise are generally lower than 1.5, also conhrm- 
ing the expected be haviour of cooling of shocked material 
dihro fc Nakar|[2013l) . 

At an effective wavelength of 4722A, the characteristic 
p-band rest-frame wavelength, we find a combined average 
rise-time duration of (7.5 ± 0.3) days. The distribution of 
rise-times spans a one a range of 4 —17 days. We caution the 
reader that given the cadences of our surveys, there could be 
more SNe II with even shorter rise-times. From this continu¬ 
ous distribution, we find no evidence for two populations of 
SNe IIP and SNe IIL. We do identify several SNe Iln and pe¬ 
culiar SN 1987A-like SNe II based solely on the long duration 
of their rise. This opens up interesting capabilities for the 
early photometric identihcation of SNe II and other SNe in 
current and future wide-held transient searches such as DES 
dBernst e in fc Dark En e rgy Survey CollaboratiorJl2014l ) and 
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Table 5. Median and deviation on the median for post-maximum slopes of the SDSS-SN and SNLS silver and golden samples at different 
effective wavelegnths. 


Effective wavelength 

SDSS-SN 

SNLS 

(A) 

silver 

golden 

silver 

golden 

3300 

3800 

4200 

4700 

5600 

6700 

8200 

U.UOZ_Q 

u.uoy_Q Q45 

U.U1D_Q QQ4 

0 m 4 + 0-006 

U.U14_q QQ4 

n m n+o-oio 

U.U1U_0.003 

U.UD4_Q 

U.U44t_Q 

U.U1Z_Q QQ3 

U.U1Z_Q QQ3 

0 m n+0-010 

U.UiU_Q QQ2 

0.043t°;°“ 

0.0l8t°;°J? 

o-oiiig'.rs 

0.008t°;°°® 

0.007t°;°J° 

0.028l°;°i4 

0.0l9l°;°i| 

0 OOQ+0-005 

u.uuy_Q QQ4 

U.UUD_q Q03 

f) QQy+0.038 

U.UU/_Q QQ4 


We compare our findings to predictions from analyti¬ 
cal models by NSIO and RWll, as well as hydrodynamical 
models from T09. The rise-times in the models are most 
sensitive to radius. The measured rise-times we obtain are 
much shorter than typical progenitor models of RSG with 
radii larger than 500 Rq. Using the analytical models and a 
new set of hydrodynamical models with smaller radii, we find 
that the bulk of our SNe are consistent with RSG progenitors 
of 200-500 7?©. The distribution of rise-times that we obtain 
translates into 84% of radii lying between 100-900 Rq. The 
mean radii we obtain are inconsistent with typical expe cta- 
tions of RSG dUevesaue et al.l[2005l : iNeilson et al.ll201ll 'l. In 
order to double the size of the radius to 8007?© in the ana¬ 
lytical models, one would need to systematically increase the 
rise-time by 7-8 days, which is well beyond our systematic 
uncertainties. These smaller radii may confir m that previ¬ 
ous R SG temperatures were under-estimated dPavies et al.l 
[201^ or perhaps th at mixing length theory is inaccurate 
dPeng fc Xioneir200ll '). 

From the shallow post-maximum slopes that we mea¬ 
sure for our sample compared to hydrodynamical models, 
we infer that large envelope masses of h ydrogen are n eeded, 
as has been previously known for SNe II dPopovlll993l i. How¬ 
ever, given the small radii found through the rise-time calcu¬ 
lations, RSG progenitors with dense envelopes are required. 
One way to achieve this is thr ough larger mixing lengt h pa¬ 
rameters in RSG progenitors dMaeder &: MeviiCTll987l 'l. We 
include new hydrodynamical models of exploding RSG with 
dense envelopes that are still not capable of explaining our 
results, but there is evidently a large range of pre-SN models 
that are not considered here. 

Alternatively, the fast rise-times and long plateaus 
could be simultaneously explained if RS G have more ex¬ 
tended envelopes or dense r oute r layers dW ittko wski et all 
I2OI2I : I Arrovo-Torres et'al] l2015l : lOhnaka et al.l l2013l l. or 
suffer mass-loss p rior to explo s ion, as has also been re¬ 
cently suggested (iGezari et al.l I2OI5I : IValenti et al.l I2OI5I : 
ISmith et al.ll201a ~). The very early rising light-curve in this 
case would not be due to shock-heated cooling but would 
be the delayed pr olonged shock break out emission of typical 
interacting SNe jMoriva et al.ll201lll . This hypothesis will 
need to be tested with corresponding theoretical models. 

Such interpretation brings up the question of the dif¬ 
ficulty to simultaenously detect the RSG shock breakout 
precursor prior to the SN rise and th e rise itself. Previ¬ 
ous claims of shock breakout detections llGezari et al. I l2008l: 
ISchawinski et al.|l2008l : iGezari et al.ll201^ were not accom¬ 
panied with a clear separation from the main light-curve. 
This could be an observational bias or could signify that such 


precursor actually does not exist because it is smeared into 
the main light-curve. Moder n high-ca dence searches such 
as Kepler (lOlling et al.ll20l3 ~l or HITS llForster et al.l[20l3 . 
I2OI5II will be able to elucidate this question with higher pre¬ 
cision than ever before. 
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APPENDIX A: NON-LTE EFFECTS ON THE 
RISE 

The numerical model of T09 was computed with the 
code Stella assuming Local Thermodynamical Equilibrium 
(LTE) conditions. Under this assumption, opacity, which 
determines the radiative transfer through a supernova en¬ 
velope, depends only on local medium properties such as 
density, chemical composition and temperature. Although 
the temperature of electrons was calculated in self-consistent 
simulations with the non-equilibrium spectrum of radiation, 
in reality the state of excitations and ionizations depends 
not only on this temperature but also directly on the de¬ 
tails of the non-equilibrium spectrum of radiation. In the 
inner envelope layers the radiation is trapped and the LTE 
assumption is fulfilled. However, at the outer edge of the en¬ 
velope, where the radiation is free to leave a supernova, the 
LTE assumption is violated. Thus the equation of state must 
be solved taking into account that the radiation spe ctrum 
is not in equilibrium with the matter (lMihalaslll978l L One 
needs to check how much the non-LTE effects will affect the 
shape of the light curve and consequently the rise time. 

To evaluate this effect, a model from the set of small 
radii by T09 was chosen. On the rising part of the light curve, 
three points were chosen at times t = 1,3,5 days (see left 
Fig. Ell • At these epochs the opacity was computed taking 
into account the non -equilibrium radiation (for details, see 
iBaklanov et al.ll2013l ). The LTE and non-LTE opacities are 
shown in right Fig. EB In the range of optical bands the non- 
LTE opacity increases comparatively to the LTE case due 
to increased Thomson scattering. This growth is associated 
with a higher level of ionization of the medium in the field 
of non-equilibrium radiation, for which color temperature is 
higher than the local electron temperature in this volume. 

Thus at the rising stage of the light curve the difference 
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Figure Al. Left: UBVRI light curves of SN models. The dotted lines show three time points t = 1,3,5 days used for the analysis of 
the role of the non-LTE effects. Right: The wavelength-dependent expansion opacities, k, are plotted at time moments t = 1,3,5 days. 
The opacity is calculated using both LTE (dotted) and non-LTE (solid) assumptions to demonstrate the reliability of the LTE approach. 


between the two approaches is small and for the integral light 
curves its effect is even smaller. In addition, the Thomson 
scattering is independent of frequency, and therefore it does 
not change the shape of the light curves. Therefore, these 
calculations show that within the accuracy of observations, 
the use of the LTE models is correct for measurements of 
the rise time. 


APPENDIX B: SYSTEMATIC EFFECTS 
B1 SNII spectral template 

We show in figure [BTla comparison between the spectral 
template by Nugent at maximum that we use in this analysis 
llNueent et al.| 2002|) wi th ear ly observed spectra of nearby 
SNe II ( Gutierrez et al.ll2014l) . Effective wavelengths calcu¬ 
lated with the observed spectra result in a wavelength shift 
of 24, 20, 64 and 2 A for ugri with respect to the template. 


B2 Cadence 

As shown in § 15.21 the SDSS-SN and SNLS have different 
cadences. This large difference leads to important system¬ 
atic biases when measuring the rise. To investigate this effect 
further we generate mock light-curves of different rise-times 
based on the model of Bazin (eq[2]) by varying the cadence 
between f and 10 days. We then calculate the rise-times for 
these light-curves in the ways mentioned in section 2] The 
rise-times based on initial dates from the power-law fits are 
generally accurate, but since we require data in the rise, 
there is a lack of objects with short rise-times when the 
cadence is high. This explains the missing fraction of SNe 
with short rise durations in the SNLS distribution. If, in¬ 
stead, we use non-detections to estimate the initial day of 
the rise, this effect is only partially mitigated (figure IB2I) 
since the end point of the rise also depends on the cadence. 
Moreover, the rise-times in this case are over-estimated, par¬ 
ticularly for short rise-times. We conlcude that this method 
does not represent an improvement over the power-law fits. 



Figure Bl. SN II spectral template by Nugent (black) compared 
to observed rest-frame MW-corrected SN II spectra of SN 2003bl 
at ~2 days after explosion (blue), SN 2008il at ~3 days after 
explosion (red) and SN 2007od at 6 days after explosion (green). 
Broadband SDSS filters are shown as dashed lines. 


The final rise-times we calculate for the SDSS-SN are always 
shorter than for the SNLS, independently of the way we 
measure them. This difference persists even at longer wave- 
lenghts, demonstrating that the distribution of rise-times in 
red bands also extends to values lower than ~ 10 days. In 
fact, the median rise-time at 8200A for the golden sample of 
the SDSS-SN is lO.llj ® days. This systematic effect opens 
up the possibility of a fraction of SNe II with even shorter 
rise-times than the lowest cadences of the SDSS-SN of 2-3 
days. 


B3 Differences in starting date 

We compare in Figure [B^ the difference between rise-times 
measured using a starting date from the power-law fit 
(eq. 0, and with a starting date from the midpoint 

of the last non-detection and first detection, 
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Figure B2. Simulated change in rise-time, measured 

with an initial date obtained from the last non-detection and first 
detection, as a function of rise-time, trise? due to varying cadence 
in the survey. The various colors indicate cadences between 1 and 
10 days. The error bars are the standard deviation divided by N, 
where = 20 is the number of simulations randomly changing 
the epoch of the first non-detection point. 

both cases, the endpoint of the rise is the maximum from 
eq. m so that all differences are due to the starting date. In 
general, we find good agreement except at short rise-times 
where the power-law measurements are larger on average by 
1-2 days for both surveys. We attribute this to the fact that 
fewer points are used in the fit, in most cases only one given 
the cadence of our surveys, and the first detection in the 
rise is closer to the time of half the flux at maximum, i.e. 
the limit of the time range used in the fit. This pushes the 
explosion date to later times. The few outliers are errors in 
the starting dates based on non-detections, for which fluc¬ 
tuations in the noise make actual non-detections sometimes 
inconsistent with zero in one particular band. This leads 
to much earlier starting dates and much longer rise-times. 
Given this systematic, we prefer using starting dates based 
on power-law fits. If the shortest measurements are slightly 
biased towards longer rise-times, the main conclusions of 
this paper are only further strengthened. 

B4 Differences in end date 

Here we compare the differences between rise-times mea¬ 
sured using the date of maximum, and using the end 

date, trSei but the same starting date. Figure [B4] shows that 
they are remarkably similar since they are both calculated 
with the same fit of Bazin (eq.[2J. However, one can see that 
for some light-curves at longer wavelengths, the maximum 
occurs during a rising plateau and is therefore quite different 
from the endpoint of the rise based on when the derivative 
becomes lower than O.Olmag/day. We thus use this date in¬ 
stead of the maximum in all Figures as a more characteristic 
end date for the rise but we quote the rise-times for the max¬ 
imum in Table [3l 

Our main findings of short rise-times leading to smaller 
radii through analytical models of NS 10 and RWll, and hy- 
drodynamical models of TOO, still hold even after putting 
the extremes of these systematics together. We cannot ac- 



Figure B3. Comparison of the effect of different starting dates 
used in the rise-time estimation: based on power-law 

fits versus ^n^det.max non-detections for the samples of 

SDSS-SN (orange circles) and SNLS (blue squares). Symbols with 
grey contours represent the golden samples and larger symbols 
are the median for both surveys. The dashed line indicate a 1:1 
relation. 



Figure B4. Comparison of two rise-times measures, i™j,^^versus 
samples of SDSS-SN (orange circles) and SNLS (blue 
squares). Symbols with contours represent the golden samples. 

count for a systematic offset of 7-8 days to recover typical 
RSG radii of ~ BOOR© in the analytical models. 


APPENDIX C: SPECTRA OF LONG RISE SNE 

Here we present the spectra of the two long rise SNe of sec¬ 
tion as a reprentative case of this set of objects. The 
observed spectra in observer-frame compared to some SNID 
templates are shown in figure ICll For SN 2006ix the better 
match to SN Hn is evident. SNLS-07D2an is an interesting 
SN that has very long plateaus in all bands which is un¬ 
expected for typical SN 1987A objects. The long rise-time 
however and the better spectral match make this a remark¬ 
able unique peculiar object. 
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Figure Cl. Observed spectra of two long rise SNe initially identified as SNe II from the post-maximum shallow slopes in arbitrary 
flux units. SN 2006ix (left) is compared to the best SN II match from SNID as well as the best SN Iln match. SNLS-07D2an (right) is 
compared to the best SN II match from SNID as well as the best SN 1987A-like match. 


This paper has been typeset from a lATp^X file prepared 
by the author. 


© RAS. MNRAS 000.[Tlf20l 

































